QUANTUM  CHAOS  AND  MESOSCOPIC  SYSTEMS 

Fiiia'  Technical  Report 


by 

Giulio  Ca,sati:  Principal  Investigator 


Julv  1997 


United  States  Army 

EUROPEAN  RESEARCH  OFFICE  OF  THE  U.S.  ARMY 
London.  E:.igla.i)d 


CONTRACT  NUMBER  N68171-96-C-9039 

Centro  di  Cidtura  Scientifica.  "A.  Volta” 
Villa  Olmo  -  22-100  Como  -  Italy 


Ai>|jrovcd  for  Public  Release;  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


form  Apcuow^d 
OMB  No.  070-i.Ql  QS 


W4V  I/O-*.  V*  lo  14*- Of****- •*  4 

- - -  4F<^«cvo.o»«  - - - 

4^^  IIM,  94  ^  .M»— r  o, 

ii^iio/o-i-aivit.  yv.,vFw..^i«A  0(; 

1.  agency  use  only  (ttjvc  C/Jn«j 

2.  RfPOBT  0Af£ 

7j~ijiiL  Amm - 

J.  REPORT  TYPE  ANO  OaTES  COVE.RfO 

Final  JReport  Januarv-March  1997 

4.  TITLE  .ANO-SUariTLE . .  .  •  -  *.  * 

Quantum  Chaos  and  Mesoscopic  Systems 

fUNCiNG  NUMBERS 

N68171~96-C-9039 

7S^2-E£-of 

■ 

6.  AUIM0R(5) 

Giulio  Casati 

Italo  Guarneri 

7.  ?£RF0KMING  ORGANIZAriCN  NAMt(S)  ANO  AOOR£SS(£S) 

Centro  di  Cultura  Scientifica  "A.  Volta" 

Villa  Olmo,  Via  Cantoni,  1  -  22100  Como,  Italy 

Centro  Sistemi  Dinamici 

Via  Lucini,  3  -  22100  Como,  Italy 

3.  performing  CPC  a. n*.;  a  tic. n 

REPORT  .VUMSER 

DYSCO  79 

I  9.  SPCNSCRInG*  .MC.VlTOPmC  AGZ.S'CY  NA.M£(S}  AnO  AOORESSitS)  j 

:o.  SPONSORING .  v.CNt: c.r;.vo  I 

w.  suppiiMcNrA^Y  notes 


;;S  r-c.v ;  a 


AfflprOTsd  ks  fgabhe  nieoKf 


:::3.  D<s :;^:au::CN  czzi 


O,  abstract  JOO^CfOz} 

The  intet-play  of  classical  chaotic  properties  and  of  quantum  coherence 
effects  is  crucially  important  in  mesoscopic  physics,  as  it  gives  rise  to 
transport  fluctuations  and  quantum  localization.  We  have  studied  both 
effects  in  different  models.  The  results  we  have  obtained  are  potentially 
relevant  to  all  fields  of  Quantum  Chaology,  with  possible  experimental 
in  the  domains  of  Atomic  Physics ,  Mesoscopic  Physics ,  and 
Classical  wave  propagation  (e.g.  acoustic,  microwave  and  optical  resonators!) 


U.  SUBJECT  TERMS 

IS.  .NUMBER  OF  paces 

21  ' 

T6.  PRICE  CODE 

IT.  SECURITY  CLASSIFICATION 

OF  report 

IS.  SECURITY  CLASSiflCATIC.N 

OF  THIS  PACE 

:9.  security  classification 

OF  ASSTAaCT 

’0.  ljmita::cn c.' abstract 

NSM  7S-0-O1.2S0-S5O0 


Stoncjro  roir\  'S3  --oS) 


ATTACHMENT  1 
PAGE  1  OF  2 


Table  of  Contents 


Abstract  . page  2 

1.  Background  . 3 

2.  Statement  of  Problems . 5 

3.  Results:  qualitative  summary  . 7 

4.  Results:  technical  description  . 8 

4.1.  Conductance  . 8 

4.2.  Billiards . 12 

5.  Conclusions  and  recommendations  . 16 

Figure  Captions  . 18 

Bibliography . 20 


1 


Abstract 


The  interplay  of  classical  chaotic  properties  and  of  quantum  coherence 
effects  is  crucially  important  in  mesoscopic  physics,  as  it  gives  rise  to  trans¬ 
port  fluctuations  and  quantum  localization.  We  have  studied  both  effects 
in  different  models.  Concerning  conductance  fluctuations  in  mesoscopic  de¬ 
vices,  we  have  exploited  a  Band  Random  Matrix  model  to  investigate  the 
lively  debated  question,  of  the  connection  between  the  different  definitions 
of  conductance,  which  are  due  to  Thouless  and  Landauer,  respectively.  We 
have  found  that  they  coincide,  apart  from  a  numerical  factor,  in  the  metallic 
regime,  but  not  in  the  insulator  one,  where  they  scale  differently  with  the 
sample  length. 

Concerning  localization,  we  have  found  evidence  that  this  phenomenon 
-  hitherto  known  for  electrons  in  disordered  solids,  or  for  externally  driven 
systems  -  is  instead  a  much  more  general  occurrence  in  the  quantal  dynamics 
of  chaotic  systems.  We  have  studied  the  classical  and  quantum  dynamics 
of  a  billiard,  which  is  chaotic  in  the  classical  limit,  and  which  is  accessible 
to  experimental  investigation.  We  have  found  evidence  for  a  quantum  lo¬ 
calization  effect,  which  sets  some  limitations  to  the  applicability  of  Random 
Matrix  theory,  and  we  have  given  quantitative  conditions  for  that. 

These  results  are  potentially  relevant  to  all  fields  of  Quantum  Chaology, 
with  possible  experimental  applications  in  the  domains  of  Atomic  Physics, 
Mesoscopic  Physics,  and  Classical  wave  propagation  (e.g,,  acoustic,  microwave 
and  optical  resonators). 
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1  Background 

Quantum  chaology  is  the  area  of  theoretical  physics  which  is  concerned  with 
the  investigation  of  the  properties  of  quantum  systems  which  are  chaotic  in 
the  classical  limit [1] [2].  It  was  started  as  a  highly  theoretical  subject,  but  has 
recently  come  closer  to  direct  physical  applications  thanks  to  developments 
in  the  field  of  Mesoscopic  Physics  [3].  This  name  describes  submicron-size 
systems  in  which  electrons  behave  quite  differently  than  in  atoms,  which 
are  much  smaller,  or  in  macroscopic  solid  samples,  which  are  much  larger. 
Typical  examples  are  conductor  or  semiconductor  devices  (e.g.,  transistors) 
of  a  few  thousand  angstroms,  which  behave  in  a  peculiar  way  due  to  quantum 
coherence  effects;  for  example,  they  exhibit  Mesoscopic  Fluctuations. 

The  fast  progress  of  microelectronics  in  manifacturing  devices  of  very 
small  size  makes  the  understanding  of  such  phenomena  a  very  important 
task  not  only  for  theoretical  physics  but  also  for  technology. 

Viewed  from  the  standpoint  of  Classical  Mechanics,  the  electronic  dy¬ 
namics  relevant  to  mesoscopic  physics  is  usually  chaotic.  On  the  other  hand, 
quantum  coherence  is  extremely  important  on  the  mesoscopic  scale,  so  that 
the  behaviour  of  mesoscopic  systems  is  actually  dominated  by  quantum  ef¬ 
fects.  The  combination  of  classical  chaotic  dynamics  and  quantum  coherence 
has  phenomenological  manifestations  in  the  form  of  localization  effects,  and 
of  fluctuations  of  transport  coefficients.  Localization  and  transport  fluctu¬ 
ations  are  therefore  a  major  feature  of  mesoscopic  physics,  and  their  study 
is  a  primary  theoretical  task,  because  they  appear  to  be  a  rather  general 
feature  of  quantum  mechanics  in  the  presence  of  classical  chaos,  and  not  just 
of  particle  motion  in  a  disordered  potential. 

Classical  chaotic  systems  exhibit  fluctuations,  which  are  not  due  to  exter¬ 
nal  random  agents  (such  as  noise,  or  thermal  effects),  but  rather  to  uncon¬ 
trollable,  small  variations  in  the  preparation  of  the  state,  which  are  rapidly 
amplified  by  the  dynamics.  On  the  other  hand,  quantum  systems  exhibit 
specific  quantum  fluctuations,  which  are  not  due  to  any  such  effect  but  are 
intrinsic  to  the  very  concept  of  a  quantum  state.  Such  quantum  fluctua¬ 
tions  are  washed  out  by  the  process  of  quantum  averaging,  which  produces 
expectation  values  that  are  usually  dynamically  stable.  There  are,  however, 
cases  in  which  even  quantum  averages  themselves  fluctuate:  such  fl,uctuations 
have  a  dynamical  origin,  and  are  presumably  connected  to  classical  chaotic 
fluctuations. 
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Phenomenologically  relevant  instances  are  well  known  from  several  fields. 
Mesoscopic  conductance  fluctuations  of  small  metallic  samples  at  low  temper- 
ature  are  a  very  popular  example[3].  Great  progress  has  been  made  possible 
in  the  study  of  such  fluctuations  by  the  use  of  Random  Matrix  Theory,  which 
approaches  the  study  of  some  relevant  properties  starting  from  the  assump¬ 
tion  that  quantum  Hamiltonians  or  Scattering  Matrices  can  be  written  as 
matrices  with  randomly  generated  elements. 

The  primary  objects  of  interest  in  the  Random  Matrix  theoretical  ap¬ 
proach  are  the  statistical  properties  of  the  energy  spectrum,  which  describe 
how  are  the  eigenvalues  distributed,  and  their  sensitivity  to  changes  in  exter¬ 
nal  parameters,  such  as,  e.g.,  external  magnetic  fields.  In  the  RMT  approach, 
a  detailed  description  of  the  spectra  is  given  up;  only  statistical  features  are 
considered.  These  turn  out  to  be  to  some  extent  universal;  in  fact,  their  gross 
features  are  already  reproduced  if  the  Hamiltonian  is  replaced  by  a  Random 
Matrix,  i.e.,  by  a  matrix  with  randomly  generated  matrix  elements,  which 
does  not  retain  any  of  the  specific  details  of  the  original  hamiltoiiian,  other 
than  its  symmetry  class  . 

According  to  current  understanding,  recourse  to  Random  Matrix  Theory 
is  justified,  when  the  system  is  quantally  ergodic.  By  this  one  qualitatively 
means,  that  eigenstates  are  maximally  extended  over  the  whole  energetically 
allowed  energy  range,  in  such  a  way  as  to  give  rise  to  phase-space  distributions 
which  fill  the  whole  classically  allowed  energy  shell.  However,  it  is  not  by 
any  means  obvious  that  such  a  structure  of  eigenstates  should  be  typical  of 
systems  which  are  strongly  chaotic  in  the  classical  limit;  in  fact  deviations 
from  ergodicity  may  occur,  due  to  several  quantum  effects. 

The  most  important  such  effect  is  dynamical  localization.  This  is  an 
interference  effect  essentially  similar  to  the  Anderson  localization  which  is 
undergone  by  quantum  particles  moving  in  a  random  potential.  Anderson 
localization  has  a  counterpart  in  the  quantum  dynamics  of  periodically  driven 
systems [4],  which  is  liable  to  experimental  detection  in  atoms  or  molecules 
in  external  electromagnetic  fields[5][6].  Dynamical  localization  results  in 
strongly  nonergodic  eigenstates,  and  severely  affects  quantum  transport,  ei¬ 
ther  in  configuration,  or  in  momentum  space,  or  in  both;  in  addition,  it  leads 
to  statistical  properties  of  energy  spectra,  quite  at  variance  with  Random 
Matrix  predictions.  , 
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Statement  of  Problems. 


As  discussed  above,  transport  fluctuations  and  localization  are  extremely 
important  issues  in  Mesoscopic  Physics,  and  their  study  was  the  main  goal 
of  the  work  summarized  in  this  report. 

2.1  Transport  Fluctuations 

There  are  two  theoretical  approaches  to  the  definition  of  the  conductance  of 
a  mesoscopic  system,  due  to  Thouless  and  Landauer  respectively.  They  will 
be  reviewed  in  the  following  section.  The  connection  between  these  different 
definitions  is  a  theoretically  intricate  problem,  as  yet  not  completely  solved. 

In  Thouless’  definition,  a  metallic  sample  is  considered  to  be  a  closed 
system,  and  the  effect  of  the  external  world  is  embodied  in  the  boundary 
conditions  to  which  the  electronic  wave  function  is  subjected.  Conductance 
is  then  related  to  the  sensitivity  of  eigenvalues  to  changes  in  the  boundary 
conditions. 

In  contrast,  in  the  Landauer  approach,  the  metallic  sample  is  an  open 
system,  coupled  to  perfect  conductors  which  carry  electronic  waves;  these 
are  partly  transmitted  and  partly  reflected  by  the  sample.  Conductance 
is  then  related  to  the  transmission  coefficient,  obtained  from  the  scattering 
matrix. 

A  special  kind  of  Random  Matrices,  characterized  by  a  band  structure 
(’’Band  Random  Matrices”,  BRM)  ,  provides  an  efficient  model  for  investi¬ 
gating  the  connection  between  the  two  definitions.  Such  matrices  have  been 
rather  poorly  investigated  until  recently,  and  are  in  a  sense  intermediate  be¬ 
tween  the  conventionally  used  Random  Matrices,  whose  essential  shortcom¬ 
ing  is  the  lack  of  a  meaningful  classical  limit,  and  the  dynamical  models  of 
Quantum  Chaos,  which  are  endowed  with  such  a  limit  but  are  more  difficult 
to  analyze. 

Band  Random  Matrices  exhibit  a  rich  variety  of  statistical  behaviours, 
depending  on  the  values  of  a  few  characteristic  parameters [7].  Different  pa¬ 
rameter  ranges  give  rise  to  several  different  statistical  regimes.  In  studying 
conductance  through  a  disordered  wire,  one  can  model  the  Hamiltonian  of 
electrons  inside  the  wire  by  a  BRM.  If  the  wire  is  considered  to  t^e  a  closed 
system,  boundary  conditions  have  to  be  given;  this  is  done  by  specifying  (yp, 
the  fixed  phase  change  prescribed  on  the  wave  function  from  one  endpoint  to 
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the  other  of  the  conducting  sample,  (p  appears  in  the  hamiltonian  matrix  as 
a  phase,  on  which  the  mati’ix  itself  is  periodically  dependent;  the  sensitivity 
of  eigenvalues  to  changes  of  this  phase  is  then  directly  related  to  conductance 
by  the  famous  Thouless  formula: 

_l  /  d^E 

A  \  d<p^ 

where  E  is  an  eigenvalue,  A  is  the  average  level  spacing  and  (.)  denotes 
average  over  disorder;  g  is  the  disorder-averaged  dimensionless  conductance. 


2.2  Localization  in  Conservative  systems 

One  of  the  main  modifications  that  quantum  mechanics  introduces  in  the 
classical  picture  of  chaotic  dynamics  is  “quantum  dynamical  localization” 
which  results  e.g.  in  the  suppression  of  chaotic  diffusive-like  process  which 
may  take  place  in  systems  under  external  periodic  perturbations.  This  phe¬ 
nomenon,  first  pointed  out  in  the  model  of  quantum  kicked  rotator  [10],  is 
now  firmly  established  and  observed  in  several  laboratory  experiments  [11]. 

How  general  is  localization  in  the  quantum  mechanics  of  chaotic  systems? 
For  conservative  Hamiltonian  systems  this  question  is  but  poorly  investi¬ 
gated.  The  situation  here  is  much  more  intriguing  :  from  one  hand,  in  a 
conservative  system,  one  may  argue  that  there  is  always  localization  due  to 
the  finite  number  of  unperturbed  basis  states  effectively  coupled  by  the  per¬ 
turbation;  on  the  other  hand  a  large  amount  of  numerical  evidence  indicates 
that  quantization  of  classically  chaotic  systems  leads  to  results  which  appear 
in  agreement  with  the  predictions  of  Random  Matrix  Theory  (RMT)[12],  and 
thus  appears  to  exclude  localization  effects,  because  the  latter  would  enforce 
deviations  from  RMT. 

In  previous  works  of  ours  [13]  the  problem  of  localization  in  conserva¬ 
tive  systems  has  been  explicitely  investigated.  In  particular,  on  the  base 
of  Wigner  band  random  matrix  model,  conditions  for  localization  were  ex¬ 
plicitely  given,  together  with  the  relation  between  localization  and  level  spac¬ 
ing  distribution  [14]. 

This  is  however  an  abstract  result,  because  the  extent  to  which  the  Wigner 
model  reproduces  the  behaviour  of  actual  quantum  systems  is  not  precisely 
known.  One  therefore  needs  to  study  more  realistic  models  of  chaotic  sys- 
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terns,  to  check  whether  the  prediction  of  localization  in  conservative  systems 
which  was  obtained  from  the  analysis  of  the  Wigner  model  is  actually  valid. 

Such  models  are  provided,  e.g.,  by  billiards,  which  are  very  important  in 
the  study  of  chaotic  conservative  dynamical  systems;  in  fact  they  are  clear, 
mathematically  well  understood  examples  of  classical  chaos,  and  their  quan¬ 
tum  properties  have  been  extensively  studied  theoretically  and  experimen¬ 
tally.  Semiconductor  devices  in  which  electrons  are  subject  to  billiard-like 
dynamics  have  been  fabricated;  moreover,  billiards  are  becoming  increasingly 
relevant  for  the  study  of  optical  processes  in  microcavities,  which  may  lead  to 
possible  applications  such  as  the  design  of  novel  microlasers  or  other  optical 
devices  [15]. 

We  have  focused  our  attention  on  a  well-known  two  dimensional  chaotic 
billiard:  the  Bunimovich  stadium,  and  have  analyzed  the  conditions  for  quan¬ 
tum  localization  and  therefore  the  conditions  under  which  the  standard  Ran¬ 
dom  Matrix  Theory  is  not  applicable. 


3  Results:  qualitative  Summary. 

Our  work  on  quantum  transport  in  mesoscopic  conductors  was  aimed  at  clar¬ 
ifying  the  lively  debated  problem  of  the  connection  between  different  defini¬ 
tions  of  conductance  of  metallic  samples  of  mesoscopic  size.  Having  numeri¬ 
cally  computed  conductance  according  to  both  definition,  we  found  that  they 
behave  differently  in  the  Ohmic  and  in  the  insulator  regime.  The  difference 
between  these  regimes  is  set  by  the  localization  effect,  and  quantitatively 
measured  by  the  localization  ratio,  i.e.  the  ratio  between  the  localization 
length  and  the  size  of  the  sample.  When  this  ratio  is  large,  localization 
effects  are  negligible,  and  conductance  scales  inversely  proportional  to  the 
sample  length,  according  to  Ohm’s  law.  In  this  regime  the  two  definitions  of 
conductance  yield  essentially  the  same  result,  apart  from  a  proportionality 
factor. 

In  contrast,  when  the  localization  ratio  is  small,  the  conductance  de¬ 
creases  exponentially  fast  when  the  sample  length  is  increased.  Here  the 
Thouless  and  the  Landauer  conductance  are  quite  different,  one  of  them 
scaling  approximately  like  the  square  of  the  other. 

The  other  question  we  have  analyzed  is,  how  general  is  localization  in 
quantum  chaotic  systems?  To  this  end  we  have  studied  the  Bunimovich 
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stadium,  a  billiard  which  is  well-known  to  be  chaotic  in  the  classical  limit 
we  have  been  able  to  show  that  quantum  localization,  until  now  believed  to 
appear  only  in  disordered  or  in  externally  driven  systems,  is  also  present  in 
conservative  systems  such  cis  billiards,  and  imposes  certain  limitations  on  the 
applicability  of  Random  Matrix  theory. 

Notice  that  the  effect  predicted  here  is  entirely  due  to  quantum  dynamical 
localization,  and  bears  no  relation  with  other,  frequently  considered,  but 
much  less  generic,  sources  of  deviations  from  RMT;  for  example,  the  existence 
of  bouncing  ball  orbits.  The  same  behaviour  will  in  fact  be  present  in  chaotic 
billiards  in  which  no  family  of  periodic  orbits  exists. 

The  effects  of  quantum  localization  discussed  here  should  be  observable  in 
microwave  or  sound  wave  experiment.  We  also  would  like  to  mention  that  the 
diffusive  process  in  angular  momentum  and  the  corresponding  suppression 
caused  by  quantum  mechanics  may  be  of  interest  for  a  new  class  of  optical 
resonators  which  have  been  recently  proposed  [15]. 


4  Results:  technical  description. 

4.1  Transport  fluctations. 

In  our  model  the  electron  dynamics  in  the  leads  is  described  by  two  semi¬ 
infinite  BRM’s  with  the  same  bandwidth  6,  with  all  matrix  elements  inside 
the  band  equal  to  1.  The  Schrodinger  equation  is 


6-|-l 

hi^i+kUi+k  =  Eui  (1) 

k=~{b+i) 

where  hij  are  the  non  zero  elements  of  the  Hamiltonian  H,  Then,  we  can 
write  H  =  Ho  +  where  Hq  is  a  BRM  of  the  same  rank  of  H  and  with  the 
same  band  width,  but  with  all  non  zero  elements  equal  to  1;  this  decompo¬ 
sition  defines  the  potential  V  of  the  scattering  problem.  Free  waves  in  the 
leads  (which  are  solutions  of  (1)  with  hij  —  1  inside  the  band)  are  of  the  form 
Um  =  with  the  wave  number  k  obeying  the  dispersion  law 


sin(l) 


^E. 


(2) 
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The  left  hand  side  of  (2)  is  a  periodic  function  of  period  27r,  which  takes 
the  value  26  +  1  when  its  argument  is  a  multiple  of  27r.  Eq.(2)  has  always 
26  solutions  for  —1<E<  +1,  which  can  be  found  analytically  for  E  =  1, 
E  =  0  and  E  =  —1  or  with  any  required  precision  by  numerical  methods  at 
any  other  value  of  E.  For  example,  they  are  kj  =  with  j  =  1  •  •  •  26  for 
E  =  0  and  ^  with  s  =  1  •  •  •  6  —  1  and  kt  =  with  t  =  0  •  •  •  6  for 
E  =  1.  The  various  allowed  values  of  the  wave- number  k  ai  a,  given  energy 
E  correspond  to  different  scattering  channels.  We  have  thus  a  multi-channel 
scattering  problem  with  exactly  6  ingoing  and  6  outgoing  channels,  which,  if 
1  «  b  «  L  ,  formally  reproduces  the  case  of  a  quasi-one  dimensional  wire. 

The  scattering  properties  of  the  conductor  are  described  by  a  unitary 
scattering  matrix  5,  that  relates  incoming  and  outgoing  amplitudes,  Il(R) 
and  Ol(R): 


where  the  subscripts  L  and  R  stand  for  the  left  and  the  right  lead.  The  S 
matrix  can  be  written  as 


S  = 


r 

t' 


(4) 


where  t  and  are  the  transmission  sub-matrices  in  the  two  opposite  direc¬ 
tions,  and  r  and  r'  are  the  reflection  sub-matrices. 

The  Landauer  conductance  is  given  by 


27rh 


El'uf. 


(5) 


where  tij  are  the  elements  of  the  t  sub-matrix.  We  have  computed  the  S- 
matrix  by  standard  methods,  by  numerically  solving  the  Lippman-Schwinger 
equation 

-  Gj  •  Vu^  =  u  (6) 

where  u  are  free  eigenfunctions,  interacting  ones,  and  the  potential  V  is 
the  difference  between  the  interacting  Hamiltonian  and  the  free  one.  The 
free  Green  function,  Gq  =  {E  —  Ho±ie)~2Q,  has  matrix  elements 


{Gt)n,m  =  {n\{E  -  Ho±ie)-^m)l^o 
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~»(m— n)fc 


E±t€- 


sin(^^fc) 


c=0 


that  were  numerically  computed  by  exploiting  the  Theorem  of  Residues. 
Finally,  the  scattering  matrix  was  computed  via  the  usual  formula: 


Sij  —  27r  \ 


Idki 

Idkj 

\l  dH 

H=E^  dH 

<  Ui \V\uf  > 


H=E 


(7) 


The  unitarity  condition  S'^S  —  1  provides  a  check  on  the  precision  of  the 
results,  and  was  fulfilled  within  a  typical  error  10“'°. 

We  have  considered  the  dimensionless  residual  conductance 


9  = 


(8) 


and  investigated  the  scaling  properties  of  the  geometric  average  Qav  —  exp  (In  g) 
in  the  regime  1  <C  6  <C  T,  where  L  is  the  length  of  the  disordered  sample. 
The  brackets  ()  mean  here  average  over  different  realizations  of  the  potential 
(1000  realizations  in  our  computations)  while  the  bar  means  average  over 
the  energy  range  (—1, 1).  As  expected,  Qav  turned  out  to  depend  on  b  and  L 
only  through  the  localization  parameter  x  —  ^  (Fig.l).  For  x  ^  1,  i.e.  in 
the  delocalised  regime,  this  dependence  has  the  form  Qav  oc  a'  oc  which 
is  the  ”ohmic”  behavior  (Fig.  2).  Instead,  in  the  localised  regime  a:  ^  1, 
Qav  -  exp(-c/a;)  (Fig.  3). 

Our  results  can  be  summarized  by  the  following  interpolating  law,  valid 
in  all  regimes: 

a  +  6exp  [—djx) 

Qav  ~  7  /  \  \ 

exp  (c/.t)  —  1 

where  c  ^  2.70  and  d  ^  0.55  (obtained  from  the  fitting  in  the  localised  re¬ 
gion)  and  a  ^  6.58  and  b  ^  —5.53  (obtained  from  the  fitting  in  the  metallic 
region).  If  the  conductance  is  computed  without  performing  the  energy  aver¬ 
age,  a  similar  scaling  law  is  obtained,  though  with  slightly  different  constants 
(curves  a  and  b  in  Fig.  4). 

We  turn  now  to  the  Thouless  approach.  In  order  to  define  the  (dimension¬ 
less)  Thouless  conductance,  the  BRM  which  describes  the  disordered  sample 
has  to  be  ’periodicized’.  The  ’periodic’  BRM  thus  obtained  describes  a  ring- 
shaped  conductor,  and  depends  on  a  phase  (f)  that  has  the  meaning  of  an 
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Aharonov-Bohm  flux  switched  on  through  the  ring.  The  disorder-averaged 
Thouless  conductance  is 


(10) 


where  A  is  the  mean  level  spacing,  and  ()  denotes  statistical  averaging.  The 
scaling  properties  of  K  were  analyzed  in  detail  in  ref.  ([7]),  which  is  the  source 
of  the  data  used  below.  In  that  work,  the  eigenvalues  sequence  was  unfolded 
to  have  unit  density,  so  that  K,  =  Kav  =  exp  (In  if),  where  K  =  |£^"(0)|.  The 
numerical  data  show  that  Thouless  conductance  obey  the  same  scaling  (9), 
but  with  constants  a'  fs  0112,  b'  ft:  — 3.12310~^,  c' m  1.565  and  d'  0.989. 

One  key  problem  in  comparing  the  Thouless  and  the  Landauer  conduc¬ 
tance  is  connected  with  the  matching  of  the  ’’free”  and  of  the  "disordered” 
part  of  the  Hamiltonian,  i.e.,  with  the  effect  of  "contacts”,  which  was  dis¬ 
cussed  in  refs. [8].  If  the  arbitrary  scale  factor  a  appearing  in  the  disordered 
part  of  our  hamiltonian  is  varied,  while  keeping  the  free  hamiltonian  un¬ 
changed,  the  Thouless  conductance  is  not  by  any  means  affected,  nor  is  the 
value  of  the  localization  length  in  the  ’infinite’  sample.  The  Landauer  con¬ 
ductance,  instead,  is  to  some  extent  modified.  It  becomes  very  small  both  at 
small  and  at  large  values  of  a,  and  it  has  a  relatively  broad  maximum  around 
a  =  1.  A  value  of  a  close  to  1  is  also  found  on  requiring  that  the  local  den¬ 
sity  of  states  for  the  ’’free”  dynamics  is  the  same  as  the  density  of  interacting 
states  in  the  sample  (having  computed  the  latter  from  the  interacting  Green 
function,  in  a  standard  way).  Our  choice  of  a  =  1,  which  maximizes  the 
transmission,  is  very  similar  to  the  "matching  wire”  prescription,  suggested 
by  Economou  and  Soukoulis. 

The  comparison  between  Landauer  and  Thouless  conductance,  computed 
as  described  above,  is  shown  in  Fig.4.  In  particular-,  in  the  metallic  regime 
(a;  >>  1)  we  have  gav  ~  (7.5  ±0.4)/'f<ji, -1-0.8,  namely  the  two  conductances  are 
proportional  in  the  diffusive  region,  where  (10)  is  assumed  to  be  valid.  The 
meaning  of  the  coefficent  of  proportionality  fti  7.5  ±  0.4  is  not  clear  to  us. 
It  is  of  course  possible  that  this  coefficient  becomes  closer  to  1  in  the  limit 
b  — r  oo,L  — r  oo,P/L  =  const.]  this  we  were  unable  to  check  due  to  obvious 
numerical  limitations,  but,  in  the  parameter  range  we  were  able  to  explore, 
the  variation  of  this  factor  is  quite  slight.  > 

In  the  localised  regime  (x  <  1)  we  get  ln(5rai;)  ~  (1-73  ±  0.09)  ln(A'at,)  + 
5.66,  that  to  say,  the  Landauer  conductance  is  proportional  to  with 
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P  =  1.73  ±  0.09.  The  large  error  in  this  numerical  estimate  is  due  to  the 
difficulty  to  compute  curvatures  in  the  localised  regime.  In  this  connection, 
we  have  also  studied  two  different  tridiagonal  models:  a  BRM  with  b  = 
1  and  an  Anderson  model  with  Gaussian  disorder  on  the  main  diagonal. 
Comparing  Landauer  conductance  and  curvature  on  these  two  models  we 
got  ^brm  —  2.078  ±  0.036  and  (S^nd  =  2.088  ±  0.0051.  It  seems  therefore 
reasonable  to  conjecture  /?  =  2;  it  is  also  worth  mentioning  that  theoretical 
arguments  suggest  that  the  dissipative  conductance  is  proportional  to  the 
squared  curvature  in  the  deeply  localised  regime[9]. 

4.2  Localization  in  billiards. 

We  consider  the  motion  of  a  particle  having  mass  m,  velocity  v  and  elastically 
bouncing  inside  the  stadium  shown  in  Fig  5.  We  denote  with  R  the  radius 
of  the  semicircles  and  with  2a  the  length  of  the  straight  segments.  The  total 
energy  is  E  =  mv^j2.  The  statistical  properties  of  the  billiard  are  controlled 
by  the  dimensionless  parameter  t  ==  afR  and,  for  any  e  >  0,  the  motion 
is  ergodic,  mixing  and  exponentially  unstable  with  Lyapunov  exponent  A 
which,  for  small  e,  is  given  by  [16]  A  ~ 

For  the  analysis  of  classical  dynamics,  a  typical  choice  of  canonical  vari¬ 
ables  is  (6,u<)  where  s  measures  the  position  along  the  boundary  of  the 
collision  point  and  Vt  is  the  tangent  velocity.  These  variables  however,  are 
quite  difficult  to  treat  from  the  quantum  point  of  view.  For  this  purpose  it  is 
convenient  instead  to  consider  /,  the  angular  momentum  calculated  with  re¬ 
spect  to  the  center  of  the  stadium,  and  the  angle  6  which  describes,  together 
with  r(0),  the  position  of  the  particle  in  the  usual  polar  coordinates.  It  is 
important  to  stress  that  with  this  choice  of  variables,  the  invariant  measure 
dfi  =  dsdvt  is  preserved  only  to  order  e,  that  is  dfi  =  dsdvt  =  dOdl  +  o(c) 

At  a  given  energy  £,  the  angular  momentum  must  satisfy  the  relation 
K|  <  Imax  =  {R  +  a)y/2Em.  It  is  therefore  convenient  to  introduce  the 
rescaled  quantity  L  =  Ill-max-  Then  the  classical  motion  takes  place  on  the 
cylinder  0  <  0  <  27r,  — 1  <  L  <  1. 

It  is  expected  that  for  e  <<  1  a  diffusive  process  will  take  place  in  angular 
momentum  with  a  diffusion  coefficient  D  =  D[t).  In  order  to  obtain  an 
estimate  for  D{e)  we  now  derive  an  explicit  expression  for  the  boundary  map 
in  the  variables  (T,  0).  The  change  AT  after  a  collision  with  the  boundary 
can  be  easily  obtained  to  order  e,  by  neglecting  collisions  with  straight  lines 
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and  by  taking  into  account  that  in  the  collision  only  the  normal  velocity 
Vn  =  V  ■  n  changes  the  sign.  Here  n  ~  +  esm9sign{cos9)eB,  eV  and  6$ 

being  the  usual  polar  coordinates  unit  vectors.  One  then  get  : 

AL  =  L  —  L  =  —2esm9sign{cos9)sign{L)\/T'—~L^  (11) 

On  the  other  side  the  change  in  9,  to  zero  order,  is  given  by  : 


A9  =  9  —  9  —  TT  —  2arcsin(//)  (12) 

According  to  a  standard  procedure  [17]  we  introduce  a  generating  function 


G{L,9)  in  such  a  way  that  the  map  defined  by 

dG 


^  89' 


9  = 


dL 


(13) 


coincides  with  AL  at  first  order  in  t  and  with  A0  at  zero-th  order.  The 
generating  function  is  given  by: 


G{L,  9)  z=z  {0  -\-  'k)L  ^  2  y  dL  arcsin  L  +  tg{L)\  cos  9\ 


(14) 


where  g{L)  =  2sz^n(L)(l— The  generated  (implicit)  area-preserving 
map  is 


L  =  L  —  2esm0sign{cos  9)sign{L){l  —  . 

6  =  6  +  TT  —  2  aTcsm{L)  +  eg\L)\ cos  6\  ^  ' 

By  taking  the  local  approximation  in  the  angular  momentum,  the  map 
(15)  writes  : 


L  ~  L  —  2esin  6sign{cos  9)sign[L)\J\  —  Lq 
0  ^  +  TT  —  2  arcsin(Z) 

which  remains  area-preserving  and  can  be  easily  iterated  (here  Lq  is  the 
initial  angular  momentum). 

The  agreement  of  map  (16)  with  the  true  dynamics  can  be  numerically 
checked  and  it  is  shown  in  Fig.6  where  we  plot  L*  =  {L  —  L)/{2eyJl  —  Lq) 
against  6,  Points  represent  billiard  dynamics  while  the  full  line  is  tlje  function 
f(^0)  =  —  sin9sign{cos9). 
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Notice  that  the  function  f{6)  is  periodic  of  period  tt  and  has  a  disconti¬ 
nuity  at  ^  =  7r/2.  This  gives  to  the  map  (16)  a  structure  very  close  to  the 
sawtooth  map  which  is  known[18]  to  be  chaotic  and  diffusive  with  a  diffusion 
rate  D  which,  for  small  values  of  the  kick  strength  e,  is  given  by  D  ^ 

This  behaviour,  according  also  to  our  numerical  computations,  appears  to 
be  generic  for  maps  which  have  such  type  of  discontinuity. 

We  proceed  now  to  a  numerical  investigation  of  the  diffusive  process.  To 
this  end  we  consider  a  distribution  of  particles  with  given  initial  Lo  and 
random  phases  0  in  the  interval  (0,27r)  and  integrate  the  classical  equations 
of  motion  inside  the  billiard.  In  Figs.  7a, b  we  present  the  behaviour  of 
AL^  =  {P)  -  (To)'  as  a  function  of  the  number  of  collisions  n  and  the 
distribution  function  fn{L)  at  fixed  n  as  a  function  of  {L  -  Lq).  As  it  is  seen, 
AL^  grows  diffusively  and  the  distribution  function  is  in  good  agreement 
with  a  Gaussian^  .  In  particular  the  dependence  D  =  D{e)  of  the  diffusion 
coefficient  can  be  easily  computed  and  the  result  D  =  (see  Fig.8)  is 

in  agreement  with  predictions  of  map  (16)  with  Do  =  1.5. 

The  analysis  of  the  classical  diffusive  process  allows  to  make  some  pre¬ 
dictions  concerning  the  quantum  motion  and  in  particular  to  estimate  the 
conditions  under  which  the  quantum  localization  phenomenon  will  take  place 
[19].  First  of  all,  in  order  that  any  quantum  diffusive  process  may  start  it  is 
necessary  to  be  above  the  perturbative  regime.  In  particular  the  level  num¬ 
ber  must  be  sufficiently  high  so  that  the  De  Broglie  wavenumber  k  of  the 
corresponding  wavefunction  must  satisfy  the  relation  k  >  1/a.  This  implies 
E  >  Ep  —  /2ma^  which  is  the  energy  necessary  to  confine  a  quantum 

particle  inside  a  box  of  length  a.  Using  the  well  known  Weyl  formula  for  the 
total  number  of  states  with  energy  less  than  E  [12] 

where  A  is  the  area  of  a  quarter  of  billiard,  and  keeping  only  the  leading 

^It  may  be  interesting  to  remark  that  the  diffusion  coefficient  computed  in  terms  of 
the  number  of  collisions,  appears  to  depend  on  the  initial  value  of  angular  momentum  L. 
This  is  due  to  the  fact  that  the  mean  free  path  depends  on  angular  momentum  and  that 
even  though  the  system  is  ergodic,  ergodicity  is  not  uniform  in  time.  If  one  computes  the 
diffusion  coefficient  in  terms  of  real  physical  time  then  the  dependence  on  tl^e  initial  Lq 
value  disappears.  For  the  above  reasons  map  (6)  approximate  the  real  dynamics  provided 
Lo  is  not  too  close  to  1. 
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term,  we  obtain  that  in  order  to  be  in  a  non  perturbative  regime  we  have  to 
consider  level  numbers 


(18) 

We  call  Np  perturbative  border. 

According  to  the  well  known  arguments  [20],  above  the  perturbative  bor¬ 
der  (18)  quantum  diffusion  in  angular  momentum  takes  place  with  a  diffusion 
coefficient  close  to  the  classical  one.  This  diffusion  proceeds  up  to  a  time 
tb  ^  Deffjh^  after  which  diffusion  will  be  suppressed  by  quantum  interfer¬ 
ence.  This  time  is  related  to  the  uncertainty  principle.  Namely,  for  times 
less  than  tb  the  discrete  spectrum  is  not  resolved  and  the  quantum  motion 
mimics  the  classical  diffusive  motion  [20,  21].  Here  D^jj  =  Doe^^^2mER^  is 
the  classical  diffusion  coefficient  in  real  (not  scaled)  angular  momentum. 

The  nature  of  the  quantum  steady  state  will  depend  crucially  on  the 
ergodicity  parameter  [21] 

V  =  ^  (19) 

te 

where  te  —  ^maxl^eff  —  2mEB? fDeff  is  the  ergodic  relaxation  time. 

For  A  <C  1  the  quantum  steady  state  is  localized  while  for  A  ^  1  we  have 
quantum  ergodicity.  The  critical  value  A  =  1  leads  to  =  Dejj  that  is 

E  =  Eerg  =  /2mR^ .  We  then  have  : 


N  =  Nerg  ^ 


1 

16^ 


(20) 


It  follows  that  only  for  N  >  N^rg  there  is  quantum  ergodicity  and  there¬ 
fore  one  expects  statistical  properties  of  eigenvalues  and  eigenfunctions  to 
be  described  by  RMT.  Instead  for  N  <  Nerg,  even  if  >  Np,  namely  very 
deep  in  quasiclassical  regions,  statistical  properties  will  depend  on  parame¬ 
ter  A  =  DoVSNe^  and  not  separately  on  e  or  N.  For  example,  the  nearest 
neighbour  levels  spacing  distribution  P{s)  will  approach  when  A  <C  1. 

We  have  tested  this  prediction  by  numerically  computing  the  level  spac¬ 
ing  distribution  for  different  values  of  e  and  N.  One  example  is  shown  in 
Fig.9  for  which  iV  >  A'p  but  since  A  <  1  the  distribution  P{s)  is  close  to 
e~‘  as  expected.  Similar  behaviour  is  expected  for  other  quantities  such  as 
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the  two  points  correlation  function,  the  probability  distribution  of  eigenfunc¬ 
tions,  etc.  The  numerical  computations  are  based  on  the  improved  plane 
wave  decomposition  method.  The  accuracy  of  eigenvalues  is  better  than  one 
percent  of  the  mean  level  spacing.  We  also  compared  the  results  with  the 
semiclassical  formula  in  order  to  check  that  there  are  no  missing  levels. 


5  Conclusions  and  Recommendations 

In  this  report  we  have  summarized  the  one-year  work  we  have  done  in  partial 
accomplishment  of  a  three  years  research  proposal.  Results  of  this  work 
appeared  in  papers[22][23]. 

In  the  study  of  electronic  transport  in  mesoscopic  conductors,  conduc¬ 
tance  is  the  primary  object  of  interest;  still,  from  the  theoretical  viewpoint, 
it  is  defined  in  two  different  ways,  which  bear  no  obvious  connection  to  each 
other,  and  we  have  supplied  empirical  evidence  for  the  relationship  between 
them. 

The  other  theme  of  our  investigation  was  localization.  Contrary  to  gen¬ 
eral  expectations,  we  have  shown  that  this  famous  effect  is  not  restricted 
to  transport  in  disordered  solids,  but  is  instead  a  rather  general  occurrence 
when  a  chaotic  system  in  which  a  classical  diffusion  occurs  is  quantized.  This 
result  has  fundamental  implications:  first,  it  indicates  certain  limits  to  which 
the  generally  believed  validity  of  Random  Matrix  Theory  as  a  description  of 
strongly  chaotic  systems  is  subject;  second,  it  brings  into  light  a  quantal 
source  of  nonergodic  behavior,  which  may  have  unexpected  implications  for 
quantum  statistical  mechanics. 

The  area  explored  in  this  research  is  potentially  relevant  to  several  fields 
of  Physics,  starting  from  those  in  which  Random  Matrix  Theory  is  recog- 
nizedly  a  basic  theoretical  tool,  and  which  include  practically  all  the  fields 
of  Quantum  Chaology,  with  direct  physical  applications  to  Atomic  Physics, 
Disordered  Solids,  and  Mesoscopic  Physics. 

Localization  effects  of  the  kind  we  have  investigated  are  important  also 
for  many  body  systems,  hence  for  Nuclear  Physics.  On  a  fundamental  level, 
they  directly  affect  the  validity  of  the  quasiclassical  approximation,  therefore 
their  study  is  crucial  in  understanding  the  quantum-classical  correspondence 
in  general,  and  the  mechanism  of  decoherence  in  particular.  Finally,  be¬ 
yond  quantum  mechanics,  they  are  also  relevant  to  classical  wave  mechanics 


16 


(electromagnetic,  elastic,  light  waves). 

The  methods  involved  in  this  investigation  range  from  standard  semiclas- 
sical  methods,  to  field-theoretical  methods  (non-linear  a-model),  to  sophis¬ 
ticated  techniques  of  numerical  simulation. 

On  the  level  of  experiments,  several  branches  of  experimental  physics 
are  involved,  including  nuclear  physics,  atomic  physics,  besides  mesoscopic 
devices. 
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Figure  Captions 


Figure  1.  Logarithm  of  the  geometric  average  of  the  Landauer  conductance 
Qav  =  vs.  logarithm  of  the  scaling  parameter  x  =  6^/L.  Nu¬ 

merical  data  (circles)  correspond  to  different  values  of  L  and  b  in  the  ranges 
25  <  L  <  400,5  <  6  <  79.  The  full  line  gives  the  fitting  law. 

Figure  2.  Geometrical  average  of  Landauer  conductance  vs.  the  scaling 
parameter  / L  in  the  delocalised  regime  »  L\  the  linear  fit  shows  the 
Ohmic  behaviour  of  Landauer  conductance. 

Figure  3.  \n{gav)  vs.  the  inverse  of  the  scaling  parameter  in  the  localised 
regime  «  L.  The  linear  fit  shows  the  exponential  decrease  of  the 
Landauer  conductance  due  to  localization:  (ln(^))  =  —CLj^^  -f  Z),  with 
C  =  2.7  ±0.04,1)  =  1.9  ±0.2. 

Figure  4.  Comparison  between  the  fittings  of  Landauer  conductance,  com¬ 
puted  with  (curve  a)  and  without  (curve  6)  energy  average  in  the  interval 
(  —  1,1).  Curve  c  shows  the  fit  of  the  Thouless  conductance  computed  in 
ref.  [19]. 

Figure  5.  The  Bunimovich  stadium  with  radius  R  and  straight  segments  2a  ; 
the  variables  {r{6)^6)  indicate  the  position  of  the  point  along  the  boundary. 

Figure  6.  Comparison  between  the  billiard  dynamics  and  the  map  (16). 
Here  we  plot  the  variable  L*  versus  6  (see  text).  Points  are  obtained  from 
numerically  integrating  the  motion  of  one  particle  in  the  billiard  for  100 
iterations,  starting  from  To  —  0  and  a  random  position  along  the  boundary, 
while  the  full  line  is  the  function  f[6)  (see  text).  Here  t  =  0.01.  The  points 
not  belonging  to  the  curve  are  due  to  collisions  with  one  of  the  straight  lines; 
this  occurrence  is  outside  the  approximation  of  map  (16). 

Figure  7.  Diffusion  in  angular  momentum  for  the  billiard  with  t  =  0.01. 
Here  an  ensemble  of  10^  particles  was  chosen  with  initial  To  —  0  and  random 
position  along  the  boundary,  a)  as  a  function  of  the  number  of  collisions 
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n;  the  dashed  line  is  the  best  fit  and  gives  D  =  AL'^jn  =  1.5  •  10"’^. 

b)  Distribution  function  after  n  =  500  collisions  averaged  over  the  last  50 
collisions.  The  full  line  is  the  best  fitting  Gaussian  with  average  —0.016  and 
variance  0.1. 

Figure  8.  Diffusion  coefficient  D  =  AL^  jn  for  the  stadium  (full  circles)  as 
a  function  of  e.  Open  circles  indicate  the  diffusion  rate  obtained  from  the 
map  (6).  The  line  is  obtained  by  the  usual  best  fitting  procedure  to  the  true 
dynamics  (full  circles)  and  gives  D  =  Doc^'^  with  Dq  =  1.5. 

Figure  9.  Level  spacings  distribution  computed  on  2000  levels  in  the  interval 
51000  <  N  <  53000  for  e  =  0.01  (a)  and  e  =  0.1  (b).  In  the  first  case 
(a)  Np  600  and  N^rg-  2.8  •  10®  and  therefore  Np  N  Nerg-  The 
value  A  ~  0.01  of  the  ergodicity  parameter  accounts  for  the  fact  that  the 
numerically  computed  P{s)  is  close  to  (full  curve).  In  the  case  (b)  one 
has  Nerg  ~  2.8  •  10®  <C  N  and  therefore,  as  expected,  the  distribution  P{s) 
is  close  to  Wigner-Dyson  (dotted  curve). 
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